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CP-PACS and JLQCD collaborations are carrying out a joint project of the 2+1 flavor full QCD 
simulation. Gauge configurations are generated for the non-perturbatively (9 (a) -improved Wilson 
quark action and the Iwasaki gauge action using PHMC algorithm at three lattice spacings, a 
0.076, 0.010 and 0.122 fm, with a fixed physical volume (2.0 fm)^. We present analysis for the 
light meson spectrum and quark masses in the continuum limit, which are determined using data 
obtained from the simulations at the two coarser lattices. Our simulations reproduce experimental 
values of meson masses. The ud and strange quark masses turn out to be m^j{ll = 2 GeV) = 
3.34(23) MeV and mf^{ii = 2 GeV) = 86.7(5.9) MeV. We also show preliminary results at our 
finest lattice spacing for which simulations are still being continued. 
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size 


a [fm] (/T-input) 


a [fm] (0 -input) 


trajectory 


1.83 


16^ X 32 


0.1222(17) 


0.1233(30) 


7000-8600 


1.90 


20^ X 40 


0.0993(19) 


0.0995(19) 


5000-9200 


2.05 


28^ X 56 


0.0758(48) 


0.0755(48) 


3000-4000 



Table 1: Simulation parameters. The production run at j8 = 2.05 is still in progress. 



1. Introduction 

The calculation of the light hadron spectrum and quark masses is the most fundamental issue in 
lattice QCD simulations. So far, systematic studies made in quenched [|l|] and two flavor {Nf = 2) 
full QCD [Q] revealed that 1) the 0(10%) deviation of the quenched spectrum from experiment 
is largely reduced in Nf = 2 QCD, and 2) dynamical up and down quarks reduce significantly 
the quark masses. As a next step in this direction, the CP-PACS and JLQCD collaborations are 
carrying out a 2+1 flavor {Nf = 2 + 1) full QCD simulation project [|3|], in which degenerate up and 
down ("light") quarks and a strange quark are treated dynamically. With this simulation, we hope 
to obtain the spectrum and quark masses much closer to those of the QCD with no approximation. 

In this article, we present the results in the continuum limit extrapolated from two lattice 
spacings together with preliminary results obtained by on-going simulations at our finest lattice 
spacing. In this study, we employ the Wilson quark formalism as in our previous works, because 
this formalism has no ambiguity in quark-flavor interpretation. We note that a similar attempt is 
being made [Q] with the staggered quark formalism. 

2. Gauge configuration generation 

For the lattice action, we employ the renormalization group (RG) improved Iwasaki gauge ac- 
tion and the clover quark action with the improvement coefficient csw determined non-perturbatively 
for the RG action [^. The choice of the gauge action is made to avoid a first-order phase transition 
(lattice artifact) observed for the plaquette gauge action 

Configurations are generated with the Polynomial Hybrid Monte Carlo (PHMC) algorithm. 
(See Ref.[Q] for our implementation.) The time step 5t in the molecular dynamics and the order 
Npoiy of Chebyshev polynomial used for an approximation of quark determinant are chosen such 
that the HMC and the global Metropolis acceptance rate achieves 85% and 90%, respectively. 

Simulations are performed at three values of the coupling constant chosen so that the square 
of the lattice spacing is placed at an even interval. The physical volume is fixed at (2.0 fm)^. 
Simulation parameters are listed in Table |l|. At each coupling, we generate configurations for 
ten combinations of hopping parameters (k",,^, Kg), five for the ud quark mass taken in the range 
mps/my ~ 0.6 — 0.78 and two for the strange quark mass chosen around mps/niv ~ 0.7. Gauge 
configuration generation has already been finished at the two coarser lattices. We plan to continue 
the on-going simulation at the finest lattice up to 6000 trajectories. 
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3. Measurement and analysis 

Measurements are performed at every 10 trajectories. We first fix tiie configuration to the 
Coulomb gauge and then calculate quark propagators for valence quark masses taken equal to one 
of sea quark masses, using point and exponentially smeared sources and sinks. For analysis, we use 
the combination of smeared source and point sink, because with this combination effective masses 
reach a plateau earliest and the statistical error is the smallest. Meson masses and quark masses are 
determined from single mass %^ fits to correlators {P{t)P{0)), {V{t)V{0)) and (A4(f)f (0)), where 
P, V and denote the pseudoscalar, the vector and the one-loop 0(a) -improved axial-vector 
current, respectively. We include correlations in time but ignore correlation among correlators. 
Errors are estimated by the binned jackknife method with a bin size of 100 trajectories. 

Chiral fits are made to light-light, light-strange and strange-strange meson masses simultane- 
ously ignoring their correlations, using a quadratic polynomial function in terms of the sea and 
valence quark masses. ^ For chiral fits, we use two definitions of the quark mass. One is the vec- 
tor Ward identity (VWI) quark mass defined by m^^^ = {1/k - \/k,)/2, where Kc is the critical 
hopping parameter at which mps at K"„^ = Kg = Kyai = K"c vanishes. The other is the axial-vector 
Ward identity (AWI) quark mass defined by m'^^' = V ^A^{x) / {2P{x)). The lattice spacing, the 
meson masses and m^^^ at the physical point are determined from the chiral fit in terms of m^^\ 
while m'^^' at the physical point are from the chiral fit in terms of m^^'. The lattice spacing and 
the meson masses obtained from the latter chiral fit are consistent with those from the former. 

For the physical point, we consider two cases. One is called as "/T-input" in which we take the 
experimental values m„ = 0.1350 GeV, ntp = 0.7684 GeV and mk = 0.4977 GeV as inputs. The 
other is "0-input" in which niji, mp and = 1.0194 GeV are taken as inputs. The lattice spacings 
determined from the K- and ^- inputs are consistent with each other, as shown in Table |l|. 

4. Light meson spectrum 

We evaluate the meson masses in the continuum limit extrapolating the data linearly in 
at two coarser lattices, for which runs and measurements are completed. The results shown in 
Fig. H are consistent with experiment though errors are still large. Taking K* with /T-input as 
an example, the linear fit is given by nin* = mo{l +c{Aqcd '(^Y) with wiq = 897(10) MeV and 
c = —1.38(78) where we assumed Agco = 200 MeV. The 0(1) magnitude of the coefficient c is 
reasonable, albeit apparently larger scaling violation is somewhat disappointing, in comparison to 
quenched(triangles) and Nf = 2(squares) data [Q] which are based on only the tadpole-improved 
action. The points at the finest lattice spacing are based on measurements for 3000 trajectories. 
Statistics are still low. We hope we continue simulations up to the full 6000 trajectories and clarify 
the continuum extrapolation in the near future. 

5. Light quark masses 

The physical quark mass is determined for the MS scheme at the scale pL =2 GeV. To do 

'chiral fit is one of the important sources of systematic errors. Recently, Wilson chiral perturbation theory (WChPT) 
has been proposed, in which infrared chiral logarithms and finite lattice spacing corrections to them for the Wilson quark 
action are incorporated. An application of WChPT to our Nj = 2+\ QCD data is in progress. (See Ref. |^.) 
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Figure 1: Continuum extrapolation of meson masses, compared with those for quenched and Nf = 2 
QCD [^. Note that the quenched and Nf = 2 simulations are made with the one-loop perturbatively 0(a)- 
improved clover action. Thus extrapolations are made linearly in a. 
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Figure 2: Continuum extrapolations of the up and down quark mass defined by the AWI and VWI (/T-input). 



this, we first translate the quark mass to that in the MS scheme at scale jj. = using Z-factors 
determined by tadpole-improved one-loop perturbation theory [^]. The renormalized quark masses 
are then evolved to jU = 2 GeV using the four-loop RG-equation. 

In principle, the quark mass calculated in our way can have 0{g^a) scaling violation since we 
use the matching factor and the improvement coefficient ca of the axial-vector current determined 
by one-loop perturbation theory. We compare in Fig. ^ the ud quark masses calculated with the 
VWI and AWI definitions. They do not agree at all in the continuum limit, when they are extrapo- 
lated in a. On the other hand, pure quadratic fits yield reasonably consistent values. This suggests 
that terms linear in a are small for ud quark masses. Similarly, the pure quadratic extrapolations of 
the VWI and AWI strange quark masses give a common value in the continuum limit. 

We take the AWI definition to quote the continuum limit for estimates of the quark mass, since 
the magnitude of scaling violation is smaller for nf^^' than for ni^^', and thus the error in the 
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Figure 3: Continuum extrapolations of the up, down and strange quark masses obtained with the A'-input. 
The data at the finest lattice is not included in the continuum extrapolations. For comparison, results for 
quenched and Nf = 2 QCD are overlaid. 



continuum limit is smaller for mf^', as shown in Fig. ^.^ Note that this property is also observed 
in Nf = 2 QCD. 

In Fig. |3[ we present continuum extrapolations of the up, down and strange quark masses 
determined with the /T-input. The ud quark mass with the 0-input agrees with the /T-input at each 
lattice spacing. For the strange quark mass, the K- and the (/(-inputs give different values at finite 
lattice spacings, but masses are extrapolated to a common value in the continuum limit, as observed 
in Nf = 2 QCD. Thus we estimate the quark masses in the continuum limit from a combined fit to 
data with the K- and the -inputs, 

m^(;U = 2 GeV) = 3.34(23) MeV, mf^{n = 2 GeV) = 86.7(5.9) MeV. (5.1) 

In the continuum limit, the quark masses estimated for the Nf = 2+1 QCD do not show deviations 
beyond statistical errors from those of Nf = 2 QCD. 



6. Conclusions and future plans 

In this article we have reported the status of our Wilson-clover Nf = 2 + 1 simulations. The 
spectrum at two lattice spacings, when extrapolated to the continuum limit are consistent with ex- 
periment. However, the outcome of the on-going runs at the finest lattice spacing at a ~ 0.076 fm, 
which we hope to complete in a half a year or so, is needed for a precise conclusion. The same 
comment also applies to light quark masses. A shortcoming of the current runs is a relatively large 
value of the dynamical up and down quark masses. We hope to bring improvement in our pro- 
gram by a combined application of the improved algorithm provided by the domain decomposition 
idea [ [To| ] and an enhanced computing power to be provided by PACS-CS, the successor to the 
CP-PACS computer. 

^The VWI quark mass for the ud quarks is negative at our simulation points. This originates from a lack of chiral 
symmetry of the Wilson quark action. This is another reason to prefer the AWI definition. 
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